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ABSTRACT 



This thesis presents an analysis of a nonlinear servo- 
mechanism. The nonlinear control device is an electro- 
mechanical spring clutch designed by the Propeller Division 
of the Curtiss-V/right Corporation for use in both missiles 

f 

and manned aircraft. 

A mathematical representation of the system is first 
derived. Then an analytic technique is presented which may 
be employed to predict such closed-loop response character- 
istics as gain, phase shift, and bandwidth. 

Investigations are made regarding both the absolute 
stability and the relative stability of the system, however, 
only limited conclusions are reached due to the lack of 
established stability criteria for a system of this type. 
Recommendations are made for further study of the stability 
of this system for the purpose of establishing such criteria. 

The ground work for this thesis, including the prelimi- 
nary analysis and all of the equipment tests, was accomplish- 
ed by the authors during the s^lmmer of I960 while on an 
industrial experience tour at the Missiles and Space Division 
of the Lockheed Aircraft Corporation, Sunnjrvale, California. 
The authors are extremely grateful for having had the oppor- 
tunity to do this work at Lockheed Missiles and Space Division 
Particularly, appreciation is expressed to Mr. Fred A. Schupan 
of Lockheed, for his very able guidance and assistance through 
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out the conduct of the analysis. Appreciation is also 
expressed to Messrs. Richard Lawhorn, Kenneth Alder, and 
Charles Cox, of Lockheed, for their assistance on various 
aspects of this work. 

The authors gratefully acknowledge the advice and assis= 
tance given by Professor Robert D. Strum of the United States 
Naval Postgraduate School during the writing of this paper. 
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CHAPTER I 



INTRODUCTION 



General 

In the field of control engineering a search is 
continually being made for new and improved methods 

f. 

for optimizing system response. The missile Industry 
in particular has presented the control engineer with a 
challenge which will tax his Ingenuity to the utmost if 
he is to supply the control devices to meet the exacting 
demands of the field. 

Various methods of control are presently in use 
throughout industry, the most common of which are a-c 
servomotors, d-c control motors, hydraulic controls, 
pneumatic controls, contactors or relays, and clutches 
and brakes. Any one of these methods may very well offer 
definite advantages over all the others in a particular 
application. Many of these devices have found wide appli- 
cation because they are linear in nature and thus are well 
suited to straightforward linear mathematical analysis. 

On the other hand, the development of other of these de- 
vices which are nonlinear in nature has been considerably 
slower because of the complexity of nonlinear mathematical 
analysis. It is pointed out, however, that due to the ex- 
tremely rapid growth of the control engineering field, and 
the common prevalence of nonlinearities in control systems. 
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the nonlinear problem has received a great deal of attention 
in recent years, and it is no longer the nemesis to the de- 
velopment of nonlinear control devices that it was a decade 
ago. 

Purpose of Paper 

It is the purpose of this paper to present an analysis 
of a nonlinear servomechanism. This nonlinear control device 
is an electromechanical spring clutch designed by the Pro- 
peller Division of the Curtiss-Wrlght Corporation. As early 
as 1952 a general study was made by Stuelpnagel and Dallas [l] 
emphasizing the advantages of a system employing basically the 
same principles that are incorporated in this system. The 
electromechanical spring clutch is, however, entirely new in 
design and construction, and Its Improved characteristics will 
serve to further advance the state of the art for clutch type 
servomechanisms . 

Scope of Study 

This paper presents an analysis of the electromechan- 
ical spring clutch functioning as the control element in a 
single-loop unity feedback system. A combination of tech- 
niques are employed in the conduct of the analysis. The 
analytic technique is used quite extensively since it is 
this approach which leads most naturally to suggestions for 
means to improve system performance through the use of com- 
pensation, More specifically, two mathematical approaches 
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are utilized. Both make use of Kochenburger > s [2] describing 
function. First, a method devised by Levinson [3] is used to 
deteimiine the closed~loop response of the system to a range 
of inputs. The Levinson technique yields system output mag- 
nitude and phase versus frequency. The system bandwidth, 

being a characteristic of primary interest in any servo system, 
« 

is readily discernible from a plot of this information. The 
other analytic approach employed is one wherein the describ- 
ing function and the linear transfer ftmction are plotted on 
a Nyquist diagram. Prom this plot limit cycle data and other 
system stability information may be obtained. 

A second technique which may be very useful in any anal- 
ysis is the trial and error approach on the analog simulation 
of the system. The analog computer is most useful in the ab- 
sence of the actual equipment itself. Since the actual clutch 
was available for this study, the analog technique found only 
limited use. 

The final technique employed in this analysis is the 
conduct of a series of tests on the system itself. These 
equipment tests, in addition to spelling out exactly what 
the system can do performance-wise, serve as a check on the 
validity of the mathematical representation used in this 
analysis . 

Some Advantages of Clutch-type Servomechanisms 

Clutches are widely used in control systems that range 
in size up to thousands of horsepower. Clutch-type servo 
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actuators compete with a~c servomotors in sizes up to one 
horsepower and with d-c servomotors over the whole power 
spectrum. One of the most Important considerations in 
servomechanism design is the speed of response. High speed 
response is difficult to obtain with a typical servomotor 
control device because of the Inertia of the motor armature. 
However, when using a servo clutch driven by a continuously 
running motor, such as the one under discussion in this 
paper, the inertia of the motor armature tends to improve 
system response time, rather than hinder it. The servo 
clutch actually utilizes the energy stored in the armature 
of the continuously running motor to give the system a high 
starting acceleration. Comparative acceleration character- 
istics of a typical servo clutch, and typical servo motor 
are shown in Figure 1. It is seen that the servo clutch 
accelerates at 30,000 radians per second per second or five 
times faster than the typical servo motor. Though not shown 
in this figure it is pointed out that the electromechanical 
spring clutch accelerates at 130,000 radians per second per 
second, or more than 20 times as fast as the servo motor. 

The electromechanical spring clutch, like most other 
servo clutches is actuated by an electromagnet. With this 
in mind, another major advantage of the servo clutch is 
pointed out. The magnetic clutch is a very effective trans- 
ducer of power from electrical to mechanical form. It trans- 
forms and amplifies power in the same operation. A comparison 
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is made between a typical 300-vxatt output servo motor and 
a 300-watt output servo clutch. In the case of the servo 
clutch, 300 watts of mechanical output power may be con- 
trolled by 10 watts of electrical input power. A 300-watt 
output from a servo motor, however, would require approx- 
imately 500 watts of controlled input power. In this com- 
pari son it is readily seen that the servo clutch alone shows 
a power amplification of 30 to 1 and an advantage over the 
servo motor of 50 to 1 in the amount of electric power re- 
quired to control 300 watts of mechanical output power. The 
electromechanical spring clutch is even more effective, for 
it requires but 2 . 21 + watts of electric input power to con- 
trol a 14.56-watt mechanical output. 
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CHAPTER 2 



DESCRIPTION OP EQUIPMENT 



General 

Before giving a detailed description of each of the 
system components, a few remarks will be made regarding 

f. 

the system as a whole. The system is an on-off modulated 
type of servomechanism. The block diagram of Figure 2 
shows the relationship of the various components. 

The signal source might be, for example, from the gui- 
dance package of a missile. The output of the servo ampli- 
fier energizes (depending upon its polarity) one of the two 
electromagnets within the clutch. One magnet causes the 
clutch to drive in one direction whereas the other causes 
it to drive in the opposite direction. Servo power is ob- 
tained from the continuously rninnlng drive motor. The drive 
motor output speed is I 4 . 23 O revolutions per minute. The cor- 
responding clutch output speed is 2^00 revolutions per min- 
ute. The linear actuator converts the rotary output motion 
of the clutch to a linear motion. Additional gearing is 
then used for angular positioning of the load, a maximtim of 
plus or minus 25 degrees from a central position. 

Photographs illustrate the actual equipment test set-up 
utilized for the analysis. See Figures 3 and Ij.. 
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Figure 2. Block diagram of on-off modulated clutch servomechanism 
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Figure 3- Laboratory set-up of closed loop system 
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Figure i;. Close-up of clutch with tachometer as used 
in laboratory tests 
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Electromechanical Sprlnpi; Clutch 

The physical size and weight of this clutch mechanism 
are Important advantages, for the clutch is but $.1 inches 
long and weighs only I |..9 pounds. 

An understanding of the electromechanical spring clutch 
is facilitated by familiarization with the nomenclature of 
the clutch assembly and examination of the sequence of events 
of operation. See Figures 5 and 6. 

The electromechanical spring clutch consists of the 
following parts: 



Clutch housing 1 
Rotary input drive 1 
Fixed shaft 2 
Rotor 2 
Electromagnet 2 



Clutch disk 2 
Clutch spring 2 
Rotary output drum 1 
Flexible shaft output 1 
Interlock mechanism 1 



Examination of the electromechanical spring clutch 
is simplified by realizing that operation clockwise or 
coxinterclockwise is identical and it is adequate to examine 
the sequence of events for a given direction of rotation. 

Use of Figure 5 and an artist's cut-away sectional drawing, 
Figure 6, greatly facilitate tracing the sequence of events. 
Mechanical input to the clutch mechanism is a constant 
speed of i|.230 revolutions per minute through gearing to the 
rotors, imparting a 2^00 revolution per minute constant speed 
to the output drum when the clutch is completely engaged and 
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Figure 5. Piinctional schematic of clutch 
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Figure 6. Cutaway view of electro-mechanical 
spring clutch 
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transients no longer exist. 

The electrical input to the clutch from the servo ampli 
fier energizes the coil of an electromagnet. The magnetic 
force attracts its associated clutch disk with sufficient 
force to hold it against the constant speed rotor. Friction 
between the clutch disk and rotor surface draws the clutch 

f, 

disk around with the rotor. The associated clutch spring is 
then tinwrapped by the clutch disk engaging a lug at the 
disk end of the clutch spring. As the clutch spring is un- 
wrapped its outer surface grabs the inner sleeve-like sur- 
faces of the driving rotor and the output drum. The clutch 
spring is maintained open with a comparatively light force 
by the clutch disk on the end of the spring. A sufficient 
torque to drive the output drum is transmitted from the driv 
ing rotor through the clutch spring. A flexible laminated 
spring shaft rigidly attached to the output drum then drives 
the output mechsinism in a rotary fashion. 

The clutch is directed to stop by removing the elec- 
trical input to the servo amplifier which de-energizes the 
coil of the electromagnet. The clutch disk then slows, 
causing the force to be removed from the lug on the end of 
the clutch spring with a resultant relaxation of the clutch 
spring. This relaxation of the clutch spring on its respec- 
tive fixed shaft causes a braking action to be Imparted 
through the clutch spring, which is rigidly attached to the 
output drum, and hence to the load. 
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The interlock gear and shaft mechanism prevents both 
ends of the clutch from engaging their respective clutch 
disks and clutch springs at the same instant and maintains 
the non=driving spring in a relaxed position. 

Linear Actuator. Gears, and Load 

The rotary output of the clutch is converted to a linear 
output through a ball-screvx type actuator. The linear output 
rate of this actuator is ^.I|.2 inches per second at the rated 
input of 2^00 revolutions per minute to the actuator. It is 
interesting to note the efficiency of this ball-screw type 
actuator is 90 percent^ compared to approximately ^3 percent 
for a comparable rotary type actuator. 

The linear output then couples directly to a rack and 
pinion type gear arrangement to which may be attached rotary 
inertia loads and friction drums. See Figure 3* 

A linear potentiometer with the pick-off tap fixed to 
the linearly driven rack gear is utilized to supply position 
feedback to the servo amplifier. 

Servo Amplifier 

Figure 7 is a block diagram showing the important fea- 
tures of the amplifier which is used to control the elec- 
trical input to the clutch. 

^Efficiency values are quoted from manufacturer's speci- 
fications. 
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Figure 7* Block diagram of servo amplifier used with 
the electro -mechanical spring clutch 



The differential amplifier produces a voltage at its 
output which is proportional to the voltage difference of the 
two inputs. It is designed to operate about a d-c voltage 
output level of seven volts. In other words, when the two 
Inputs add algebraically to zero, the output of the differ- 
ential amplifier is seven volts. The amplifier gain may be 
adjusted without varying the bias level of the output. 

Preliminary system tests were conducted using only the 
differential amplifier to excite the clutch coils. For these 
tests, the output bias level was adjusted to zero. This re- 
sulted in merely applying a proportional signal to the clutch 
colls. This type of input was found to be unsatisfactory. At 
a particular magnitude of the proportional signal, the clutch 
was foTUid to be at a point where it was about to engage, but 
the signal was not of sufficient strength to cause full or 
positive engagement Immediately. This resulted in slippage 
within the clutch which in turn caused excessive heating of the 
clutch. To overcome this problem an electronic relay, the 
Schmitt trigger, was introduced into the amplifier circuit. 

As can be seen in Figure 7 two such triggers are used. The 
operation of the two Schmitt triggers can best be understood 
by referring to Figure 8. It is readily seen that there is 
present in the firing cycle of the two triggers a slight hys- 
teresis effect. As long as it is kept small it will not se- 
riously affect the performance of the system. Hysteresis 
will be discussed more fully in a later section. The object 
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Figure 8. Schmitt trigger firing cycle 



fMli^ 




i ' ii 

i \ * 










I 








f 






I 



I 




of using the Schmitt trigger, then, is to delay the signal 
to the clutch coils until it reaches a threshold value, and 
upon reaching this value, a voltage will be applied to the 
clutch coils which is of sufficient strength to insure posi“ 
tive engagement of the clutch, thus eliminating any possi- 
bility of slippage. 

f, 

Initial open- loop tests conducted on the system indicat- 
ed the delay in stopping the clutch to be considerably longer 
than the delay in starting. Since these delays have great ef- 
fect on the response characteristics of the system it was de- 
sirable to reduce them as much as possible, particularly the 
long stopping delay. Toward this end, a reverse kick genera- 
tor was built into the amplifier. Its operation is explained 
as follows. Assume that coil number one has been energized 
for some time, and is about to be cut-off. When trigger num- 
ber one cuts off, this fires reverse kick generator number 
two, which in turn generates a signal of sufficient magnitude 
to fire Schmitt trigger number two. Clutch coil number two 
is thus energized, the associated electromagnet clamping the 
number two clutch disk, thereby creating a drag which assists 
in bringing the clutch to a stop. Obviously the duration of 
the reverse-kick signal is very important. If it is too long 
the clutch will not remain stopped but will drive in the op- 
posite direction. Thus provision was made so that the duration 
of the reverse-kick signal may be adjusted to any length be- 
tween two and 25 milliseconds. 

19 
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The power gates serve merely to supply the clutch coils 
with a signal of desired magnitude, once the Schmitt triggers 
have fired. 
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CHAPTER 111 



PERFORMANCE CRITERIA FOR THE SYSTEM 

Performance specifications in servo literature follow 
no \mlversal pattern and are as numerous as there are authors. 
Even the more completely developed linear servo field lacks 
uniformity, therefore, prior to listing the specifications it 
seems appropriate to discuss some of the major items which 
must be considered in the establishment of criteria for a sys- 
tem of this type. 

Since there is an exact correspondence between the opera- 
tion of an on-off clutch control circuit and the circuits of 
on-off relay control, [4], the criteria to be discussed below 
are precisely that prescribed in the literature, [5], for 
contactors and relays. 

In an on-off control system one of the major design con- 
siderations is static accuracy. By this it is meant, the 
range of error which can exist when the servo is in a state 
of rest without initiating corrective action. This require- 
ment is met by keeping the range of the Inactive zone within 
specified limits. The narrower the inactive zone, the greater 
will be the static accuracy of the system. However, narrow- 
ing the inactive zone tends to increase system gain and there- 
by decrease system stability. Thus it is seen that a com- 
promise between static accuracy and stability considerations 
may be necessary. 



21 







mm «»<# INV»«» 

mm- « 



4iMt 



■Ml* III 
••»! JUm^ 



A second consideration is that of dynamic accuracy, or 
the extent to which errors are minimized when the system is 
responding to a disturbance. The accuracy of a system under 
dynamic conditions is dependent upon the speed of response 
and the degree of stability. The speed of response is in 
turn, dependent upon the runaway correction magnitude. (The 
runaway correction magnitude, R , is defined as the steady- 
state rate of change of the output that would be approached 
if the clutch were to remain engaged for a sufficiently long 
period.) The higher this magnitude, the faster the system 
responds to a given disturbance. But, as in the case of de- 
creasing the inactive zone, increasing the riuiaway correction 
magnitude corresponds to an increase in gain and therefore 
generally tends to make the system less stable. Here again, 
a compromise is necessary between speed of response and sta- 
bility considerations. 

Another major consideration with regard to establishing 
performance specifications is that of system stability. What 
is the absolute stability of the system, and what is the de- 
gree of stability? As to the absolute stability, an on-off 
control system may come under one of the following categories 

1) Be completely stable. Here there is no tendency 
for self-sustained oscillations or for constantly increasing 
errors following a disturbance. 

2) Have self -sustained oscillations of a definite am- 
plitude following the application of a sufficiently large 
disturbance . 
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3) Be conditionally stable. In this case, either the 
system has no output oscillations, or such oscillations are 
limited to a given amplitude provided the disturbances to 
it do not exceed some critical value. Disturbances above 
this value can cause ^instable operation to result. 

Ij.) Be ^stable. In this case, any disturbance, no 
matter how small, results either in a constantly increasing 
error or in oscillations which increase in anqslitude toward 
destructive breakdown. 

Now, the degree of stability takes on a different mean- 
ing depending upon the category of absolute stability. For 
systems which are always stable, "degree of stability" refers 
to the rate at which the oscillations resulting from a given 
disturbance are damped. This is not to be confused with the 
"damping ratio" of linear systems. In the on-off control 
system the pattern of such damped oscillations is dependent 
upon the magnitude of the disturbance, and constants such as 
"damping ratio" and "logrlthmic decrement" have no signifi- 
cance. 

For systems which maintain self- sustained oscillations, 
"degree of stability" takes on a slightly different meaning. 
Suppose a large disturbance is applied to a system of this 
type. The initial oscillations may be of large amplitude, 
but these oscillations diminish in amplitude until the steady- 
state amplitude of the sustained-oscillation condition is 
reached. Specifications can be established in terms of the 
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time required for oscillations to be reduced to a specified 
amplitude following a given disturbance. 

Conditionally stable systems might be acceptable where 
it can be assured that no disturbance will ever occur of such 
a magnitude as to cause instability. In such a case, degree 

f- 

of stability specifications may be applied similar to those 
outlined above. 

Still another factor which warrants some consideration 
and may be prescribed in the performance specifications is 
the ntimber of engagements the clutch makes. The reason for 
this is obvious. The greater the number of engagements, the 
greater the wear on the clutch and the greater will be the 
heat generation within the clutch. 

As mentioned previously the size of the Inactive zone 
plays an important part in determining static accuracy, dy- 
namic accuracy, and stability in general. It is also ex- 
tremely important in regard to number of clutch reversals 
required to bring the clutch to a stop with zero error. The 
clutch must be stopped before it crosses the dead zone, or 
it will drive full speed in the opposite direction and will 
never be brought to a stop. Thus it is evident that both 
the Inactive zone and runaway correction magnitude must be 
studied jointly in an effort to solve the stopping problem. 

In summary, let it be said the above discussion by no 
means covers all the details that must be considered in es- 
tablishing the performance criteria for an on-off type servo. 
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It is only meant to emphasize some of the major considera- 
tions involved in such a process. It may be said then, that 
the performance requirements for the on-off system may be 
met by specifying a range of inactive zone consistent with 
static accuracy requirements and a runaway velocity consis- 
tent with thf desired speed of response. These, of course, 
must be chosen with an eye toward adequate stability which 
may or may not require some compensation effort. 

Having in mind the factors which must be considered in 
arriving at a set of specifications the following performanc 
criteria^ are listed? 

Runaway correction magnitude, R = / sec 

Width of the inactive zone, A = 1° to 2^ (range) 

These specifications define an adequate* servo system 
and are kept quite general remembering that the purpose of 
this paper is to present an analysis of this system. Other 
characteristics of the system operating under these speci- 
fications are to be determined from the analysis. Detailed 
compensation and optimization are not objectives of this 
investigation and will be mentioned only briefly in ensuing 
sections. 



^Design characteristics of the equipment test set-up, 
inactive zone, and runaway correction magnitude were 
established by Lockheed Aircraft Corporation as typical 
requirements suitable for flight control application. 



CHAPTER IV 



PRELIMINARY ANALYSIS 



General 

A series of preliminary tests were conducted on the 
open=loop system in an effort to determine accurately any 
time delay characteristics the clutch might possess, and 
to establish one other very Important point, the effect of 
load on the response characteristics of the clutch. The 
effect of load would have considerable bearing on the scope 
of the closed=loop analysis. Refer to Figure 9 for a 
diagram of the system on which these preliminary tests were 
made. The tests were conducted for various conditions of 
externally applied load at various magnitudes of command 
signal. Three magnitudes of a step function input were 
used? 28, 45, 60 volts, as measured at the coils of the 

clutch. The load conditions were 0, 20, and 45 ft-lbs of 
torque applied externally by means of a friction druun. In 
addition to this applied torque, an inertia disk was added 
at the output, resulting in another 6.2 slugs mass at the 
output of the actuator. 

Delay Characteristics 

Table I is a summary of the results obtained in the 
preliminary tests. Figure 10 shows oscillogram traces of 
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Figure 9. Block diagram of system for preliminary analysis 



Table I 



Results of Preliminary Tests 





Coil Voltage 


Load 


Condition 


Delay 










Start 


Stop 


(1) 


28 volts 


0 


ft-lb 


12 msec 


I 4.5 msec 


ftf 


20 




n 


n 




!t 


1+5 




ft 


II 


(2) 




0 




8 


1+5 




TT 


20 




ff 


It 




It 


1+5 




11 


If 


(3) 


60 


0 




8 


1+5 




If 


20 




It 


It 




If 


1+5 




11 


tf 



the position and velocity outputs for (2) above. Figure 10(a) 
shows the delay time for starting to be approximately 8 msec, 
Sind Figure 10(b) verifies the stopping time to be that given 
in condition (2) of Table I. 

These tests were made using a simple variable gain ampli- 
fier, Observing the stopping time to be undesirably long, the 
amplifier circuitry was modified to include the reverse kick 
generator in an effort to reduce this time. See Figure 7. 
Through experimentation it was found that a reverse kick of 
about 5 msec duration would give a stopping time of approxi- 
mately the same duration as the starting delay. A reverse 
kick of any longer duration would cause full reversal of the 
clutch. Then, by using a reverse kick applied to the opposite 
coil, the starting and stopping delays were made equal. This 
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Figure 10. Delay characteristics of clutch \inder various conditions 
of load and volts to coil 






will, of course, simplify the mathematical analysis of the 
closed=loop system, 

A limit cycle analysis indicated that seven milliseconds 
was a more representative value for the magnitude of the delay 
of the system, thus it was believed that the magnitude of the 
delay was reasonably well established. In the mathematical 
analysis to follow, the delay will be represented as a trans- 
portation lag, where t = 7 msec. 

Effects of Load 

The manufacturer's specifications indicate that the 

clutch is designed to respond to load torques of 25.1 in- lbs 

at the clutch output shaft without slippage. This corresponds 

to a load torque of l80 ft-lbs at the output of the system 

/ 

under analysis. Though the clutch is capable of responding to 
loads of such magnitude, the range of interest in this anal- 
ysis extends only to approximately 50 ft-lbs thus loads beyond 
this magnitude vrere not applied in the preliminary tests. Now 
the effect of load on the response characteristics of the 
clutch was determined from actual tests to be negligible. 

These results were as follows; 



Load Condition 

0 ft-lb 
20 



Acceleration Time 

1.0 msec 

1.0 

1.0 



These data may be read directly from the velocity curves of 
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Figure 10(a). Similarly no appreciable effect is observed 
in the stopping characteristics shown in Figure 10(b). 

These results indicate that for loads within the range 
of interest, the response time of the clutch is not affected 
by load. With this Information in mind, it was concluded 
that no further investigations into the effects of load on 
response time need be made. 
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CHAPTER V 



MATHEMATICAL REPRESENTATION OF CHARACTERISTICS 
Clutch Characteristics 

It was learned in the preliminary analysis that the 
clutch operates essentially as an on-off control device with 
a constant time delay. That is to say, the clutch, in re- 
sponding to a disturbance, remains at a standstill during the 
delay, after which time it is driving at full speed. It was 
established that the acceleration time in attaining full speed 
from a standstill is approximately 1.0 msec. This time is 
negligible when compared with the delay, for example. Thus, 
for analysis the clutch will be considered a nonlinear on-off 
device possessing hysteresis. 

Mathematically then, the clutch vjlll be represented by 
both a transfer function, , and a describing function, G^ . 
The transfer function will actually represent the dynamics of 
the clutch and actuator. The describing fxmction will serve 
as a mathematical representation of the on-off characteristics 
of the system. The delay previously established will take the 
form of a transportation lag, in the mathematical repre- 

sentation of the system. Figure 11 shows, in block diagram 
form, the system arrived at for the mathematical analysis. It 
will be noted that no amplifier is shovai in this diagram. The 
amplifier gain is, of course, taken into account in the magni- 
tude of the inactive zone of the nonlinearity. A decrease in 
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Figure 11, Mathematical representation of system 



the size of the inactive zone corresponds to an increase in 
amplifier gain. 

Transfer Pxmction, 

The transfer function which represents the linear com- 
ponents of the system will now be developed. As shown in 
Figure 11 it has the general form 

R(TiS + 1) 

^ s(T 2S + 1) (T^s + 1) 

where R ~ 175°/sec is the rxinaway correction magnitude set 
forth in the specifications. The T terms, which will be 
derived below, represent the clutch dynamics as well as 
those of the actuator and the rack and pinion. The linear 
factor s , in the denominator of the transfer ftinction, re- 
presents a single order pole at the origin of the complex 
s-plsine and it Identifies a system having a sustained rate 
for a constant input signal. 

To give the reader a better understanding of what is 
represented by the transfer function, a block diagram showing 
the mechanical features Involved is presented in Figure 12. 
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Figure 12, Block diagram showing mechanical featvires of system 
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The values of the various constants shown in Figure 12 are 
as follows? 

n^ = l/l „69 (drive gear ratio) 

kg = 5000 in-lb/rad 

J = 59 x 10 "^ in- lb-sec^ 
s 

k - 261 in-lb/rad 

f = 1.3 in-lb/rad/sec 

J = 78x10“^ in-lb-sec^ 
a 

n = 1/66 (actuator gear ratio) 
a 

Jq = 517.7x10"^ in-lb-sec^ 
n^ = 1 (output gear ratio) 

It will now be shown that any contribution to the trans- 
fer function that would be made by the drive spring constant, 
kg , and the inertia term, Jg , may be neglected. The fre- 
quency at which this terra might begin to contribute to the 
system response would be at 

kg 5000 

CO = — = = 9200 rad/sec . 

^s 59x10“^ 

The anticipated frequency range of interest will extend only 
up to 300 radians per second, maximum. The frequency, 00 _ , 

o 

Is so far above this range that its contribution to the sys- 
tem response will be negligible. Therefore, the terra in- 
volving k and J will not be included in the derivation of 
the transfer function. 
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The equations of motion for the remaining elements shown 
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be written 


as follows: 


(1) 
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f(01 - 0p) 
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J = J + j' 
a o 

Solving equation (i;) for initial conditions equal zero; 
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Evaluating the coefficients and factoring, the transfer func 
tion becomes 



§1 



^205 ® ^8555 



s + 1) 



By including the runaway velocity, R , and the s in the 
denominator, the total transfer fTinctlon, Gg , then becomes 



175 ( 



201 



s + 1) 



1 1 
s(tt3 s + 1) (T nm j s + 1) 



205 



8555 



Examining the significance of the various terms of this trans- 
fer function in the complex s-plane it is seen that the effect 
of the pole at s = -205 is nullified by a zero at s = -201 . 
And furthermore, since the frequency range of Interest is re- 
stricted to values below 300 radians per second, the term in- 
volving the frequency of 8555 radians per second may be dis- 
regarded, The transfer function to be used in this analysis 
then is 



G 



s 



175 

s 



Describing Function, G^j 

The describing function was developed by Kochenburger 
as a means for analyzing a contactor servomechanism. Since 
the electromechanical spring clutch functions in the same 
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on-off manner as a three-position contactor, the describing 
function will be employed in the frequency response analysis 
of this paper. 

The describing function, by definition, is the ratio of 
the amplitude of the fiindamental component of the output of 
the nonlineapity, to the amplitude of the sinusoidal input to 
the nonlinearity. Having previously determined the character- 
istics of the clutch, they are shown schematically in Figure 13 
for the purpose of establishing the notation which will be used 
in the development of the describing function. A signifies 
the width of the Inactive zone, and h , the width of the hys- 
teresis loop. The assumed input to the nonlinearity is the 
sinusoidal waveform of Figure ll 4 .(a). The resultant waveform of 
the correction signal is shown in Figure llj.(b) . Referring to 
these two figures, the describing ftinction is now derived as 
follows 2 
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1 

2 



cos 



A -h 
2 lEl 



cos 



A +h 
2 lEl 



By means of a Fourier analysis, the magnitude of the fvindamen- 



tal component of the correction signal, 
be 

i N-, i 



N, 



may be shown to 
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Figure 13 . On-off characteristics of clutch 



Control Signal, E 




(a) Assumed input to nonlinearity 



Correction Signal, N 



Actual signal 
Fundamental component 




(b) Output of nonlinearity and its 
fundamental component 



Figure 1[|.. Relationship of correction signal to control 
signal for on-off system possessing both an 
inactive zone and hysteresis 



kl 



and the phase angle, 0 , Is seen in Figure lit ( b) to be 
0 = ™ cx , 

the minus sign indicating lag. The describing function, 
is now written as 





k 

tt!e! 



sin /3 /— 



For the case of no hysteresis, the phase angle is zero, and /3 
simplifies to 

= cos~^ — ^ — 

2 |E| 

giving the rather simplified form of the describing function, 
This completes the mathematical representation of the 



system. The transfer function, the describing function and 
the transportation lag are all that is needed to proceed with 
the mathematical analysis of the closed- loop system. 



CHAPTER VI 



CLOSED-LOOP RESPONSE OP THE SYSTEM 



Analytic Development 

The closed-loop frequency response of a nonlinear system 
often is not easy to determine. Difficulties arise because 
for a given system input and frequency the amplitude of the 
actuating signal, E , depends on , but is in turn 

determined by E , As a result of this interdependence of 
G^ and E , the calculation of the closed-loop frequency re- 
sponse for a given input is at times a rather complicated pro- 
cedure. In solving the problem this analysis makes use of an 
approach which is attributed to E. Levinson [3]. A graphical 
technique is used in determining values of E which satisfy 
both frequency and nonlinear conditions. The details are ex- 
plained in the development which follows. 

Referring to Figure 11, an equation for the error signal, 
E is derived 
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Recalling 



G. 



TTlEl 



sin /3 /- OC 



let 



nlEi 



sin 



/S = 



N 



then 



also recalling 



Gj = N cos OC - jN sin OC 
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Since it is the frequency response which is of Interest, 
substitution s = jco is made, giving 

175 



G 






Similarly for the transportation lag, 
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“ cos (^2 - j sin 0j^ 

Making the substitutions equation (2) now becomes 
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The magnitude of the error may then be written as 
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|ei| O) 
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j^[l75N cos((X+ 0 ^)]^ + [co- 175N sin(cx + 0 ^)] 
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Squaring both sides of equation (5), 
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©il^ 
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'l7^N cos(cX + + [co - 175N sin(0C + 0j^)J 
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1e|2 [i75n cos(OC+ 0-^)Y 

+ |E|2 [00=^ 175N sin(oc + 0 ^)Y 
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transposing terms, 



|e!2 [i75n cos(cc + 

|2 [co- 175N sln(OC+ 



0ii . IE 



(7) 



Taking the square root of equation (7) , 



|e| [175N cqs(OC + 

= 1^ \ 9±\^U)^ - |E|2 [oj - 175 N sin(a + 



( 8 ) 



Rearranging terms and factoring, equation (8) may be 
expressed in the form 



(9) 



NE 



CO 



175 cos( OC + 0^^ L 



|9i|2 - |E|2[i - sln(OC+ 



Novj, since the input to the system is assumed to be sinusoidal 



0^ = A sin Got 



therefore, [©il^ = 



Making this substitution, equation (9) becomes 



(10) 



NE 



CO 



175 cos(OC + 0 



— Ta2 - |e|2[^i - sin( 0 C + 0^)] 

1^ L 



Each side of equation (10) is evaluated separately. A 
width of inactive zone, A , is first chosen. For the example 



whose results are plotted in Figure l5, a value of A = 2.01° 
was selected. Preliminary equipment tests indicated an essen- 
tially constant value of the ratio h/A = 0.091 . Thus, for 
the selected value of A , h = 0.183°. The left side of equa- 
tion (10) is now evaluated, letting |e 1 vary over the antici- 
pated range ,of error signal. The results are plotted in 
Figure l5 as the [NEI curve. 

To evaluate the right-hand side of equation (10), the 
value of A remains ■unchanged. A magnitude of input to the 
system is next chosen. For the exsimple whose results are 
plotted in Figure 15, A = 10.7°. There still remain two 
variables, co and |E|, in this part of the equation. The 
next step is to select a value of co , then letting |E| vary 
over a sufficient range, a curve is obtained which is also 
plotted in Figure 1^ and labeled with its associated value 
of CO, say co = 2TT . Another value of co is selected and 
the procedure is repeated yielding another curve, say the 
Co = 4TT curve of Figure l5. This procedure is repeated un- 
til the complete frequency range of interest has been covered. 

The points of intersection of the |NEl curve with the 
CO j' curves yield values of |E| which satisfy equation (10) 
and thus satisfy both frequency and nonlinear conditions. 
Having determined values of N and 1 e| over the desired fre- 
quency range the magnitude of the system output is known from 
the relationship 
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( 12 ) 



Since it is coiranon practice to express the closed-loop 
response in .form of gain and phase versus frequency, the 
gain equation is obtained by dividing equation (12) by the 
magnitude of the input, 
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Similarly, the phase of — (joj) may be obtained. Look- 

®i 

Ing at the equation for this ratio 
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175N /-OC-01 
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[175N cos( oc + 0]^)] + J [co - 175N sin( OC + 0]^)] 



From equation (l5) the phase of 
by 0 is 



!f 

0. 



■(jco) 



49 



which is denoted 



= = 



[Cij - 175 N sin( (X + !Zf^)] 

[l?5N cos(oc + !2f^)] 



(16) 



OC+ 0^ + tan 



-1 



With the aid of Figure l5, equation (16) may be evaluated 
to give the desired phase versus frequency plot. 

«, 

Comparison of Analytic and Equipment Results 

The gain curve obtained from evaluating equation (13) 
for the assumed values of 0^ and A is plotted in Fig- 
ure 16(a). For ease of comparison, on the same plot is shown 
the corresponding curve obtained from an equipment run util- 
izing the same parameters. It will be noted that the analytic 
and equipment gain curves compare quite favorably over the 
entire frequency range. 

The phase results which were obtained from the evalua- 
tion of equation (16) are plotted in Figure 16(b). Here the 
analytic approach gives extremely close correlation with 
equipment results up to a frequency of approximately 100 ra- 
dians per second, beyond which the two curves begin to di- 
verge slightly. Nevertheless, the two curves compare reason- 
ably well over the entire frequency range. 

For an input, = 10.7°sin uj t and an inactive zone, 

A = 2o01°, Figure 16(a) indicates the gain response of the 
system to be essentially flat until velocity saturation occurs 
at about 20 radians per second. The system bandwidth for the 
given conditions is approximately five cycles per second. 



50 



I 




f 




0 




Response to Other Amplitudes of Input 



The closed=loop response of the system to other magni- 
tudes of input, A , were calculated by the above method for 
comparison with the results obtained on the equipment. The 
inactive zone, A , was maintained at 2.01°. While the j NE | 
curve of Figure 15 is still applicable so long as A remains 
tmchanged, the curves which Intersect with the |NE| curve must 
be recalculated each time a new magnitude of is chosen. 

This means re-evaluating the right-hand side of equation (10), 
for each new 0£ . This was done for magnitudes of 0^ above 
and below the first example given, nsimely for A = 21.1|.° and 
A = 5.35°. The gain and phase results for 9j[^ = 21.1;° sin cot 
are shown in Figures 17(a) and (b) respectively. The same 
close correlation is seen to exist between analytic and actual 
equipment results. The bandwidth for the higher magnitude of 
input has decreased as one might expect, to approximately 2.5 
cycles per second. The response curves for an input of 5.35° 
are presented in Figures l8(a) and (b). The bandwidth in this 
case is approximately nine cycles per second. 

Response without Reverse Kick Generator and with A = t;.9° 

The gain and phase loci for the system operating without 
the reverse kick generator and with an inactive zone of 1;.9° 
are shown in Figures 19(a) and (b). To emphasize the true 
significance of the reverse kick generator in the amplifier 
circuit, a comparison is made between the curves of Figures 
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16(b) and 19(b), Figure 19(b) Is a plot of phase versus fre- 
quency of the system with the reverse kick generator removed 
from the circuit and using an inactive zone, A= l4-.9°. Fig- 
ure 16(b) on the other hand shows phase versus frequency for 
the system using the reverse kick generator, but with A= 2.01°. 
Now the effect of the wider inactive zone may be seen by com- 
paring the two analytic curves. It is seen the wider inac- 
tive zone results in a greater phase shift at lower frequen- 
cies but as the frequency increases the two analytic curves 
of Figures 16(b) and 19(b) differ less and less. Practically 
speaking the phase shift is essentially the same for the two 
magnitudes of A for frequencies above approximately four 
cycles per second. At 20 cycles per second, for example, both 
analytic curves indicate a phase shift of approximately 135 *^ 
lag. Theoretically then, the phase shift observed from the 
equipment runs for the two magnitudes of Inactive zone should 
be essentially the same beyond four cycles per second. Com- 
paring the solid curves of Figures 16(b) and 19(b) it is seen, 
however, that the system using no kick-back (Figure 19b) is 
lagging considerably behind that where kick-back is employed, 
(Figure l6b) , At four cycles per second the system without 
kick-back has a phase shift of -^9° as compared with -i|6° for 
the system using kick-back. And at 20 cycles per second the 
system without kick-back has a phase shift of -168° versus 
- 150 ° for the system using kick-back. This is certainly suf- 
ficient evidence of the importance of the reverse kick genera- 
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tor in improving the system response. 

The response gain curve for the system using the wider 
inactive zone and no kick-back is shown in Figure 19(a) . it 
Is Interesting to note, however, that the gain or magnitude 
results shown In this figure are essentially the same as 
those obtained for the system using kick-back and an inactive 
zone, A - 2,01° (Figure l6a) . 

Typical Response of the System to a Sinusoidal Input 

In the previous sections response gain and phase 
characteristics have been discussed at some length. To sup- 
plement this discussion it seems appropriate to show some 
representative examples of the system response to an applied 
sinusoidal input signal. In Figures 20(a), (b), (c), and (d) 
are shown the response to an input, 0^^ = 10.7°sln cot for 
frequencies from one cycle per second to six cycles per sec- 
ond. In Figures 21(a), (b), (c), and (d) are shown similar 
curves for an input, 0^ = 21.ti.°sin cot. The width of the 
inactive zone is 2.01° for both sets of curves. 
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(a) Frequency of one cycle per second 





Figure 20, Typical response curves for a sinusoidal 
input of amplitude 10,7^ ( A = 2.01°) 
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(d) Frequency of six cycles per second 



Figure 20, Typical response curves for a sinusoidal 
input of amplitude 10.7° (A = 2.01°) 
Continued 
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(a) Frequency of one-half cycle per second 




Figure 21, Typical response curves for a sinusoidal 
input of amplitude 21,l\P (a = 2.01°) 
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(d) Frequency of four cycles per second 



Figure 21. Typical response curves for a sinusoidal 
input of amplitude 21 . 14 .° (A = 2.01°) 
Continued 
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CHAPTER VII 



SYSTEM STABILITY 



Absolute Stability 

To examine the stability of the system, the response 
ratio is expressed as 

9o,. , 

— (JUJ) = ; 

®i 1 + G,G 

d s ^ 



G -1^ j + Q 

s ^ d 

The test for stability depends only upon the denominator 
of equation (l?) since the numerator is independent of fre- 
quency and merely represents a scale factor [ 2 ]. Thus the 
study of stability involves an investigation of the poles 
of + G,^ , or the values of frequency and ampli- 

tude which satisfy the equation 

(j^-l^jOjt = (;L 8 ) 

Substituting for the inverse transfer function, G_“^ , 

s 

equation (I 8 ) becomes 



.]■ 00 C 
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= -G. 



(19) 
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Now the procedure is to plot each side of equation (19) 
separately for the purpose of analyzing system stability. 

In the preliminary analysis the value of t in equation (19) 
was determined to be seven milliseconds. Using this value, 
the left=side of equation (19), which shall be referred to 
as the frequency locus, is plotted in Figure 22 for values 
of to up to 300 radians per second. For purposes of illus- 
tration, the value of A = 2.01° which was used in the 
closed-loop analysis of chapter VI is assumed in evaluating 
the right-hand side of equation (19). The latter, which 
shall be referred to as the amplitude locus, is also shown 
on the Nyquist plot of Figure 22. The points along the 
amplitude locus represent the magnitude of the error signal 
which is really the input to the nonlinear clutch. The re- 
lative orientation of the amplitude and frequency loci deter- 
mines the stability characteristics of the system. Since the 
amplitude locus always appears to the left of an observer as 
he traverses the frequency locus in the direction of increas- 
ing frequency the system will be stable for all amplitudes of 
input. Suppose now, the amplifier gain is increased until 
the corresponding width of inactive zone is 0.962°. The cor- 
responding amplitude and frequency loci are plotted in Fig- 
ure 23. The frequency locus, of course, is identical to that 
of Figure 21, however the amplitude locus is changed until it 
now intersects the frequency locus at two points. System 
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stability is then interpreted as follows [6j. The system 
is unstable for 0.57° < |E| < 0.85° and stable for all other 
magnitudes of error signal. If |E1 = 0.57° sustained oscil- 
lations will result. If |E| increases slightly above 0.57°, 
the oscillations will increase until sustained oscillations 
of magnitude |E| = 0.85° result. The frequency of oscillation 
is then determined by the point of intersection of the fre- 
quency and amplitude loci where |E| = 0.85° . In Figure 23 
this frequency is shown to be approximately 220 radians per 
second. A still larger increase in the control signal above 
!E| = 0.85° results in a stable system, with a decay of os- 
cillation amplitude back to |E| = 0,85°. The point on the 
amplitude locus corresponding to |E| = 0,57° represents un- 
stable equilibritim conditions and the point where |E| = 0.85° 
represents stable equilibrium conditions. It may then be 
concluded that the system \inder study, when operating with an 
Inactive zone, A = 0.962° , will oscillate continuously with 
an |e| = 0.85° and at a frequency of approximately 220 radians 
per second, provided the magnitude of the control signal, 

|E| ^ 0.57°. For 1E| < 0.57° the system is stable with no 
oscillations resulting. 

Comparison of Predicted and Observed Results 

For comparison with the predicted results outlined in 
the previous section the system as shown in Figure 11 was 
taken with zero input. This corresponds to the system 
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operating open=loop„ The amplifier gain was then increased 
xontil an inactive zone, A = 0.962° was obtained. The sus“ 
tained oscillation which resulted had an amplitude of ± 1.5l° 
and a frequency of 225 radians per second. The predicted and 
observed frequencies of oscillation differ by approximately 
2.3 percent. On the other hand, though the predicted and ob- 
served amplitudes are of the same order of magnitude, they 
differ considerably. Figure 23 shows the predicted amplitude 
to be 0.85°. 

In many applications the discrepancy that was foimd to 
exist between predicted and observed amplitudes would be of 
little significance, however in some particular application 
closer correlation might be desirable. No attempt will be 
made in this study to fully justify this difference, however 
it is pointed out that the answer may be related to the shape 
of the waveform existing at the frequency of oscillation. 
Under the stated test conditions, the output waveform is also 
the input to the nonlinear device. At these higher frequen- 
cies, this waveform is triangular rather them sinusoidal as 
ass\3med for the application of describing function theory. 

Degree of Stability 

The absolute stability of the system has been discussed 
at some length in the previous sections. Of equal importemce 
is the relative stability of the system. The system must 
possess sufficient damping so as to eliminate any undesirable 
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oscillation during the correction of a disturbance. For 
linear servomechanisms, Brown and Campbell [7 ] have estab- 
lished a rule of thumb procedure which gives an indication 
of the degree of stability from frequency response loci. A 
peak value Mp of the output-input ratio, ®o^®i deter- 

mined, The rule states that in general there is adequate 
damping if Mp < 1.3 . The rule is used widely in design 
work, however it naturally does not apply in all cases. Its 
use must be backed up by the designer's experience. 

Kochenburger developed a similar rule for use in anal- 
yzing the relative stability of a contactor servomechanism 
[ 2 I. Since the system under study functions similarly to a 
three-position contactor, Kochenburger’ s technique will be 
employed here in an effort to give the reader some feel for 
the relative stability of this system. 

In this nonlinear system there is not just a single 
value of Mp over the frequency range. There is an Mp 
for each magnitude of control signal, E , These values of 
Mp and the associated frequencies, oj p are determined 
graphically as follows. Referring to equation (I 7 ) it will 




( j '-w ) 



Gr 






(17) 



be noted that G-^ is a fixed quantity for a given magnitude 
of the control signal, therefore to determine Mp , the peak 
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value of the ratio ©q/®! ’ necessary only to determine 

the denominator of equation ( 17 ) such that it is a minimum. 
Taking as an example the system with a A = 2.01°, Figure 22 , 
which represents this system, is expanded and reproduced as 
Figure 2I4. so that the graphical work will be more legible. 

The quantity of equation (1?) is represented in Figure 2i|. 

by the vector DO , The denominator of equation (1?) is re- 
presented by the vector stun (DO + OG) or the vector DG . A 
value of M is then determined graphically from the ratio of 
vectors DO and DG , 



DG 



This, however, is the value of M corresponding to the fre- 
quency CO = 50 radians per second. M , being the peak 
value of M over the frequency range, is determined by the 
minimum value of DG . The minimum length of DG is fotind 
by drawing an arc of shortest possible radius from the point 
D tangent to the frequency locus, shown in Figure 

24 for a magnitude of |E| = 1.2°. Then merely by evaluating 
the relationship 



DO 

M == 

P DG 



min 



over the anticipated range of the error signal the loci, M 
and CO versus !E| may be obtained. These loci, for the 
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system where = 2,01°, are plotted in Figure 2^. 



Interpretation of Relative Stability 

Figure 25 indicates that the system possesses excellent 
damping characteristics at relatively low natural frequencies 
during the beginning of a correction process involving high 
amplitude ’^control signals. As the correction process nears 
completion the degree of damping decreases as revealed by a 
rise of the Mp locus, however Mp is never greater than 1.6 . 
At the lower values of |E| the frequency associated with 
transient oscillations would Increase considerably as Indicated 
in Figure 25 o 

The authors, with their limited experience, will not 
attempt to interpret the full significance of the relative 
stability of this system within the scope of this study. It 
is not suggested that one evaluate the degree of stability of 
this system according to rules established by Kochenburger for 
contactor servomechanisms, nevertheless it is interesting to 
note these criteria, and then note the results shovm in Fig- 
ure 25 for this system, Kochenburger found from experimental 

tests on various types of contactor systems that an M of 

P 

2,0 near the cut-off point and of 1.3 for higher amplitudes 
of error signal provides a satisfactory degree of relative 
stability for most applications. Figure 25 shows this system 
to have an M < 1,6 near the cut-off point and M - 1 at 
higher values of !e| , In other words, the system which 
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is the subject of this study possesses a greater degree of 
relative stability throughout the range of the control sig- 
nal than Is considered adequate for most contactor servo- 
mechanisms. 



«; 
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CHAPTER VIll 



CONCLUSIONS 



This study demonstrates quite conclusively a mathemat- 
ical technique for determining the closed-loop response char- 
acteristics of a servo system employing the Curtiss-Wright 
electromechanical spring clutch as the control device. It 
has been shown that the closed-loop gain and phase loci 
arrived at analytically by this technique agree very closely 
with the actual results obtained from equipment tests. The 
system bandwidth, which of course in a nonlinear system is a 
ftinction of the input signal amplitude, may be read directly 
from the gain versus frequency locus. 

While the calculations involved in arriving at the 
closed- loop response of a system by the method employed in 
this paper are quite straightforward, they are nevertheless, 
extremely tedious. Therefore, where extensive use of this 
technique is anticipated it is recommended that a digital 
computer be used to perform the calculations. 

The investigation into absolute system stability yielded 
close agreement between theory and actual tests in the fre- 
quency of the limit cycle oscillations. However, extremely 
close agreement was not obtained between the theoretical and 
actual amplitudes of these oscillations. The shape of the 
control signal at these high frequencies is triangular and 
it is the opinion of the authors that this may account for a 
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major part of the discrepancy in the amplitudes since the 
theory assumes a sinusoidal control signal. 

Based on the relative stability criteria for contactor 
servomechanisms, a system employing the electromechanical 
spring clutch has adequate stability. Even though both sys- 
tems function as on-off control systems, the validity of 
using such criteria is questionable, and it is believed that 
this is an area in which additional study might be done to 
definitely establish these stability criteria for a system 
employing the electromechanical spring clutch. 
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